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We review the confinement scenario in Coulomb gauge. We show that when thin center vortex configurations are gauge 
transformed to Coulomb gauge, they lie on the common boundary of the fundamental modular region and the Gribov region. 
This unifies elements of the Gribov confinement scenario in Coulomb gauge and the center-vortex confinement scenario. We 
report on recent numerical studies which support both of these scenarios. 



1. Introduction 

In a confining theory such as QCD it is helpful to 
choose a gauge which makes the confining mecha- 
nism transparent. In the present article we shall be 
primarily concerned with the confinement scenarios 
in minimal Coulomb gauge and in maximal center 
gauge, and with the unification of elements of these 
two scenarios. 

The Landau gauge has the simplest Lorentz trans- 
formation properties, but there is at present no con- 
finement scenario in Landau gauge. The obstacle is 
that the gluon propagator is of shorter range in Lan- 
dau gauge than it is in a free theory, because of the 
suppression of the low-momentum components by the 
proximity of the Gribov horizon (see below). This 
makes the mechanism of confinement of color charge 
in Landau gauge more rather than less mysterious. We 
turn instead to the Coulomb gauge. 

Confinement of color charge is easily understood in 
minimal Coulomb gauge because the 0-0 component 
of the gluon propagator, 

DQQ{x,t) = 14oui(|a;|) S{t) + non-instantaneous, (1) 
has an instantaneous part, K;oui(''), that is long range 
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and confining and couples universally to all color- 
charge. We call Vcoui('') "the color-Coulomb po- 
tential". Moreover the 3-dimensionally transverse 
would-be physical components of the gluon propaga- 
tor, 

D,,{x,t) = (A,(x,t)Aj(0,0)}, (2) 

are short range, corresponding to the absence of glu- 
ons from the physical spectrum. We shall review re- 
cent numerical evidence [1-3] which strongly sup- 
ports these statements. 

The Coulomb gauge has also been studied vigor- 
ously in the hamiltonian formalism [4]. 

2. Definition of minimal Coulomb gauge 

In minimal Coulomb gauge, the representative 
Ai{x) of each gauge orbit (at a fixed time t) is the ab- 
solute minimum of a minimizing functional with re- 
spect to gauge transformations. The minimizing func- 
tional is taken to be 

FA{9)^\\'A\f (3) 
where the gauge transform is given by 
'~'A,^g-'A,g + g-'d,g, (4) 
and the Hilbert norm by, 

/ d'xj2\An^)\'- (5) 
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The set of representatives, A, satisfies 

A^{^,(x):||A||<||M||}, 



(6) 



for all gauge transformations g{x), and is called the 
"fundamental modular region". 

In lattice gauge theory, the minimal Coulomb 
gauge is defined analogously. A maximizing func- 
tional is defined by 



(7) 



Gauge orbit \ 

Fundamental modular region , 
Gribov region 



' Along 






Atr 





and the fundamental modular region by 

A = {[/ : Fu{l) > Fuig)}. 



(8) 



The maximization is done independently within each 
time slice t. Euclidean and Minkowski Coulomb 
gauges are the same. Only physical excitations are 
propagated in the time direction. 

In practice it is difficult to find the absolute mini- 
mum, and one is satisfied to find a relative minimum. 
In principle one should check the sensitivity to choice 
of minimum. 

3. Elementary properties of minimal Coulomb 
gauge 

At a relative or absolute minimum, the minimizing 
functional is stationary, which yields the transversal- 
ity or Coulomb-gauge condition 



^,A^{x) = 0. 



(9) 



In addition, the matrix of second derivatives is non- 
negative. In the present case it is the Faddeev-Popov 
operator, 

MiA) = ~d,D,{A) > 0, (10) 

where Df^iA) = 8^5°-" + f'''"'A\{x) is the gauge- 
covariant derivative. These two properties hold at any 
relative or absolute minimum. Together they define 
the Gribov region 



VL^ {A: properties (9) and (10) hold}. 



(11) 



The set of relative minima fi is larger than the set of 
absolute minima A, and we have the inclusion A C O, 
as illustrated in Fig. 1. In continuum gauge theory 
both A and VL are convex and bounded in every direc- 
tion [5], and in SU(2) lattice gauge theory a slightly 
weaker convexity property is established in [3]. 



Figure 1 . A typical gauge orbit through the configura- 
tion Ai is represented as a curve. Transverse configu- 
rations lie in the horizontal plane, viewed from above. 
It contains Q, and A, and the pair of Gribov copies Ai 
and A2 ■ 

4. No confinement without Coulomb confinement 

The gauge-invariant, physical potential V{R) be- 
tween a pair of external quarks is obtained from a 
rectangular Wilson loop, of dimension R x T, with 
T large, 

exp[-V{R)T] = (Pexp J A^dx''). (12) 

Whereas V{R) involves n-point functions of all or- 
ders, the color-Coulomb potential T4oui(^), defined 
in (1), is the instantaneous part of Z3oo, a 2-point func- 
tion. 

The ground state energy lies lower than the 
Coulomb energy, which leads to the interesting in- 
equality [6] 



E,, + V{R)<E'^^ + V,oni{R), 



(13) 



where Ege and E'^^ are self-energies that are finite in 
the presence of a lattice cut-off. This bound is ob- 
tained from a trial wave-functional, and relies on the 
fact that the same kernel K'^y {A) whose expectation- 
value is the color-Coulomb potential, 

V,oui{\x-y\)5ab^{Kll{A)), (14) 

also appears in the Hamiltonian in Coulomb gauge for 
a static quark pair at x and y. 



H 



{{n*'^,Kfy{A)t''^) 



.(15) 
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This kernel has the expHcit form 

KfyiA) = [M-\A){-V^)M-\A)r,l, (16) 

where M{A) is the Faddeev-Popov operator (10). 
It follows that if V{R) is confining, 

lim V{R) = 00, (17) 

ii— »CX3 

then the color-Coulomb potential, Koui(^)> is con- 
fining, 

lim Vcoui(i?) = 00. (18) 

R—>-oo 

Moreover if both increase linearly at large ii, V{R) ~ 
a R; ycoui(-R) ~ CTcoui R, then 

CTcoul > (19) 

We conclude that in the confining phase the 2-point 
function Vcoui(^) is confining. 

Greensite and Olejnfk [1] determined the color- 
Coulomb potential Koui {R) numerically from the lat- 
tice quantity 

Gi\x\) ^ ^{Tv[U^{x,t)Uo{0,t)]), 
VoouiiR) = -lnG{R). (20) 

They found an impressively Unear behavior of the 

color-Coulomb potential K'oui(^)- The relevant fig- 
ure is presented in Greensite' s talk at this meeting [7]. 

5. Conimement scenario in Coulomb gauge 

To see why Vcou\{R) is long range, consider for- 
mula (16) for the kernel K^y{A). Recall that M{A) is 
strictly positive in the interior of the Gribov region fl, 
and develops a zero eigenvalue on its boundary (the 
Gribov horizon), 

M{A) > for A inside fl, 
M{A)(t>o = OforAondn. (21) 

By continuity, M{A) has a small eigenvalue for con- 
figurations near the boundary. As explained below, 
the population is concentrated close to the boundary 
dn. This enhances M-^{A), and thus also K^yiA), 
which makes Koui(-R) long range. A more sophisti- 
cated discussion, given below, explains the long range 
of V^oui {R) in the infinite- volume limit in terms of the 
density of states p{X) of the Faddeev-Popov operator 
M{A). 



As already noted by Gribov [8], the Gribov horizon 
is close by in directions (in A-space) that correspond 
to low-momentum Fourier components of the gluon 
field A. This suppresses the low-momentum parts of 
all Lorentz components of the gluon propagator D^^ 
in minimal Landau gauge. While this eliminates glu- 
ons from the physical spectrum, it makes it harder 
to explain confinement of colored quarks. For the 
same reason, in minimal Coulomb gauge, restriction 
to the interior of the Gribov horizon suppresses the 
low-momentum components of the would-be physi- 
cal, 3-dimensional gluon propagator {x, t), which 
eliminates physical gluons from the physical spec- 
trum. However, as we have just seen, in Coulomb 
gauge it also enhances V"coui(|a;|), which is the instan- 
taneous part of Dqo in Coulomb gauge. This couples 
universally to all colored objects, and is a prime can- 
didate to trigger confinement. Thus we can have our 
cake and eat it: physical gluons are suppressed be- 
cause Dij is short range, while 2?oo is long range and 
can cause confinement. 

6. Double boundary dominance 

The dimension D of A-space in the presence of the 
lattice cut-off is a large number, of the order of the 
volume of the lattice D ^ , and in a space of high 
dimension the volume density goes like r^^^dr. We 
thus expect on entropy grounds that (a) the population 
is concentrated close to the boundary OA. However, 
as we have just seen, the enhancement of Koui(-R) 
occurs if (b) the population is concentrated close to 
the boundary dil. However A is a subset of fi, A C fi, 
so (a) and (b) are compatible only if (c) the population 
is concentrated where the 2 boundaries coincide, i.e. 

dAndn. 

Why should this be? We present an explanation in 
terms of dominance of center vortices, and thereby 
unify elements of the Gribov confinement scenario 
in Coulomb gauge and the center- vortex confinement 
scenario. The latter is reviewed in [9]. 

7. Thin center vortex configurations lie on the 
double boundary 

Definition: we call a "thin center vortex configu- 
ration" any lattice configuration for which every link 
variable is a center element, Ux,i = Z^^i. The uni- 
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fication of elements of the Gribov and center-vortex 
scenarios follows from the fact [2] that when a cen- 
ter configuration Z is gauge transformed to the min- 
imal Coulomb gauge Z— >[/ = sZeAit Ues on 
the common boundary, U e dQ. n 9A. The proof is 
sketched below. It follows that center vortex domi- 
nance implies dominance by a subset of this common 
boundary. 

Note: The same statement and conclusion hold also 
for abelian configurations. 

8. Thin center vortex configurations lie on singu- 
lar gauge orbits 

Thin center vortex configurations and their gauge 
transforms U = ^Z may be characterized in a gauge- 
covariant way as possessing non-zero solutions ^„ 

D{U)(l)n = 0, (22) 

where n = 1, . . . - Ifoi SU{N), and D{U) is 
the generator of a gauge transformation, defined by 

3U = U + eD{U)(f), (23) 

where g = 1 + ecf) is an infinitesimal gauge trans- 
formation. Indeed, center configurations are invari- 
ant, = Z, under the N'^ — 1 linearly independent 
global (x-independent) gauge transformations h. For 
h = 1 + e(j) this gives the gauge-covariant condition 
D{Z)(t) = 0. 

Note: Similar properties hold for abelian configu- 
rations, with n= 1, . . . ,R, where R is the rank of the 
gauge group. 

Proof that a thin center vortex configuration Ues 
on the double boundary: The equation D(Z)<pn = 0, 
being gauge covariant, holds after transformation to 
minimal Coulomb gauge, U — ^Z, so U G A, and 
with ipn = ^(pn, we have D{U)ipn = 0. It follows 
that the tpn are nuU vectors of the Faddeev-Popov op- 
erator 

-ViA(f/)V'n = 0, (24) 

so U = ^ Z lies on the boundary of the Gribov region, 
i7 e Sil. We also have U € A, and since A c fl, 
it follows that U lies on the double boundary, U € 
dA n dn, as asserted. 

The gauge orbit of a thin center configuration is a 
geometrically singular gauge orbit. It has N'^ — 1, 
fewer dimensions than a generic gauge orbit, because 



9/7 = U for g{t) = cxp(f(^„). By contrast, the Gri- 
bov horizon dfl is, in general, merely a coordinate 
singularity. 

9. Numerical tests of center vortex dominance in 
Coulomb gauge 

In [3] the hypothesis of center vortex dominance 
was tested by the following procedure [10]. (1) Con- 
figurations are fixed to the maximal center gauge. 

(2) In this gauge, a thin center-vortex configuration 
is defined, for SU(2), by 

Z^,^ = sign[Tr([/^,^)], (25) 
and a "vortex-removed" configuration by 

UxfiJ, = Zx,ij,Ux,ij,- (26) 

(3) Finally, the resulting center-projected configu- 
ration Zj. ^ and the vortex-removed configuration 
Ux,fj. are gauge transformed to the minimal Coulomb 
gauge. According to the center vortex scenario, the 
projected center configurations should retain the con- 
fining properties, whereas the vortex-removed config- 
urations should be non-confining. 

The color-Coulomb potential K;oui(^)> defined in 
(20) was determined numerically in [1], for both the 
full configuration and the vortex removed configura- 
tion. The relevant figure is reproduced in Jeff Green- 
site's talk at this meeting [7]. One sees that the 
color-Coulomb potential is impressively linear for the 
full configuration, as noted above, but is flat for the 
vortex-removed configurations. This confirms the 
center- vortex scenario. 

A more detailed test of the center-vortex scenario 
comes from a study [3] of the density p{X) of the 
eigenvalues A of the Faddeev-Popov operator (10). 
This quantity appears in the formula for the Coulomb 
or unscreened self-energy of a static color-charge at 
infinite lattice volume that follows from (14) and (16), 

^ = ^'"°'^(p(A)F(A)), (27) 

where F(X) = (A|(— V^)|A) is the diagonal matrix 
element of (— V^) in the Faddeev-Popov eigenstates 
I A). The Coulomb energy is infrared divergent, as it 
should be in the confined phase, if, at infinite volume, 

limi(p(A)F(A))>0. (28) 
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Eigenvalue density, p=2.1 
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Figure 2. The F-P eigenvalue density at /3 
8* - 20^ lattice volumes. 



2.1, on Figure 4. F-P eigenvalue density in vortex-only con- 
figurations. 



m=(<fx i-^'i-Pa). m-1 



4 





full configurations 




^1— L=8 ^ ■ 




••^ L=10 — - 




L=12 : 




L=16 • 




L=20 



0.2 0.4 0.6 0.8 1 



Figure 3. F{X), the diagonal matrix element of 
(— V^) in F-P eigenstates, plotted vs. F-P eigenvalue. 

The quantities p(A) and F{X) were determined nu- 
merically [3] for the full configurations, for center- 
projected (vortex-only) configurations and vortex- 
removed configurations, each of which have been 
transformed to Coulomb gauge. 

Figures 2 and 3 show the results for (p(A)) and 
{F{X)) for the full configurations, on a variety of lat- 
tice volumes ranging from 8^ to 20^. The apparent 
sharp "bend" in p(X) near A = becomes increas- 
ingly sharp, and happens ever nearer A = 0, as the 
lattice volume increase. The impression these graphs 
convey is that in the limit of infinite volume, both 
p{X) and F{X) go to positive constants as A — » 0. 
However, for both p{X) and F{X) we cannot exclude 
the possibility that the curves behave Uke A^, A^ near 
A = 0, with p, q small powers. 

Next we consider the same observables for the 
"vortex-only" configurations, consisting of thin cen- 
ter vortex configurations transformed to Coulomb 
gauge. The data for (p(A)) and {F(X)) at the same 
range (8^ — 20^) of lattice volumes are displayed in 



FW=<<Pl l-V"l<Pl>,|3=2.1 
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Figure 5. F{X), the diagonal matrix element of 
(— V^) in F-P eigenstates, for vortex-only configura- 
tions. 

Figs. 4 and 5. The same qualitative features seen 
for the full configurations, e.g. the sharp bend in the 
eigenvalue density near A = 0, becoming sharper 
with increasing volume, are present in the vortex-only 
data as well, and if anything are more pronounced. 

Finally, we consider the same observables for 
the vortex-removed configurations, transformed to 
Coulomb gauge. Results are shown in Fig. 6 for 
(p(A)). The behavior is striidngly different, in the 
vortex-removed configurations, from what is seen in 
the full and vortex-only configurations. The graph of 
(/9(A)), at each lattice volume, shows a set of distinct 
peaks, while the data for {F{X)) (not shown) is or- 
ganized into bands, with a slight gap between each 
band. Inspection shows that eigenvalue interval as- 
sociated with each band in {F{X)) precisely matches 
the eigenvalue interval of one of the peaks in {p{X)). 

In order to understand these features, consider the 
eigenvalue density of the Faddeev-Popov operator 
^xy = appropriate to an abeUan the- 
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ory, or a non-abelian theory at zeroth order in the 
coupling. At finite lattice volume, this operator has 
degenerate eigenvalues, and we call the degener- 
acy of its k-th eigenvalue Afe. When the degenera- 
cies Nk, of the zeroth-order Faddeev-Popov eigen- 
values are compared with the number of eigenval- 
ues per lattice configuration found inside the fc-th 
"peak" of (p(A)), and fc-th "band" of {F{X)), there 
is a precise match. This leads to a simple interpre- 
tation: the vortex-removed configuration C/^ can be 
treated as a small perturbation of the zero-field Umit 
Ufj, = I. This perturbation hfts the degeneracy of 
the Afc, spreading the degenerate eigenvalues into the 
peaks of finite width in (/9(A)) seen in Fig. 6. For 
the vortex-removed configurations, both (/9(A)) and 
{F{X)) seem to be only a perturbation of the corre- 
sponding zero-field results. 

We conclude that it is the vortex content of the ther- 
malized configurations which is responsible for the 
enhancement of both p(A) and F{\) near A = 0, lead- 
ing to an infrared-divergent Coulomb self-energy. 

10. Conclusion 

We have seen that when thin center vortex con- 
figurations are gauge transformed to the minimal 
Coulomb gauge they lie on the double boundary of 
the Gribov region and the fundamental modular re- 
gion. This unifies elements of the Gribov confinement 
scenario in minimal Coulomb gauge and the center- 
vortex confinement scenario. 

Our numerical study [3] reveals the following fea- 
tures: 

(1) The data are consistent with a linearly rising 
color-Coulomb potential, K:oui(^) ~ o'coui R, and a 
Coulomb string tension that is larger than the physical 
string tension, ctcouI > c. 

(2) The data are compatible with a density of eigen- 
values /9(A) of the Faddeev-Popov operator with ei- 
ther p{0) finite, or with p{X) ~ A^', for small A, 
where p is close to zero. This feature, and the value 
of (p(A)F(A)) at small A, provide detailed verifica- 
tion of the Gribov confinement scenario in Coulomb 
gauge. 

(3) These confining features are preserved by the 
"vortex-only" configurations, but are replaced by fea- 
tures close to a free theory in the "vortex-removed" 
configurations. This is consistent with center vor- 
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Figure 6. F-P eigenvalue densities for vortex-removed 
configurations, on a 20* lattice volume. 

tex dominance. This, in turn, implies the condition 
of double boundary dominance that accords with the 
Gribov scenario in minimal Coulomb gauge. 

Finally, we refer to [3] for a similar numerical in- 
vestigation of the gauge field coupled to a Higgs field, 
and of pure gauge theory at finite temperature in the 
deconfined phase. 

REFERENCES 

1. J. Greensite and S. Olejnik, Phys. Rev. D 67, 
094503 (2003), arXiv: hep-lat/0302018. 

2. J. Greensite, S. Olejnik, and D. Zwanziger, 
Phys. Rev. D 69, 074506 (2004), arXiv: hep- 
lat/0401003. 

3. J. Greensite, S. Olejnik, and D. Zwanziger, arXiv: 
hep-lat/0407032. 

4. A. Szczepaniak, Phys. Rev. D 69, 07403 1 (2004), 
arXiv: hep-ph/0306030; C. Feuchter and H. Rein- 
hardt, arXiv: hep-th/0402106; D. Zwanziger, 
arXiv: hep-ph/03 12254. 

5. M. Semenov-Tyan-Shanskii and V. Franke, Zap. 
Nauch. Sem. Leningrad. Otdeleniya Matematich- 
eskogo Instituta im. V. A. Steklova, AN SSSR, 
vol. 120, p. 159, 1982 (Enghsh translation: New 
York, Plenum Press 1986). 

6. D. Zwanziger, Phys. Rev. Lett. 90, 102001 
(2003), arXiv: hep-lat/0209105. 

7. J. Greensite, These Proceedings, arXiv: hep- 
lat/0407022. 

8. V. Gribov, Nucl. Phys. B139, 1 (1978); D. Zwan- 
ziger, Nucl. Phys. B518, 237 (1998). 

9. J. Greensite, Prog. Part. Nucl. Phys. 51, 1 (2003). 

10. Ph. de Forcrand and M. D'EUa, Phys. Rev. Lett. 
82, 4582 (1999), arXiv: hep-lat/9901020. 



